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Abstract 

There are four IL-1 Receptor-Associated Kinase (IRAK) 
family members in mammals, one of which is IRAKI. IRAKI 
interacts with IRAK4, which is associated with MyD88 
adaptor protein upon ligand stimulation. This research aims 
to study the domains of IRAKI, understanding of which is 
important to yield insights into the functional role of IRAKI 
in the pathogen sensing mechanism. The results obtained 
from this study may guide further research in the 
immunotherapy studies on the Toll-Like Receptor (TLR) 
signaling pathway. 
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Introduction 

IL-1 Receptor-Associated Kinase 1 (IRAKI) is an 
important kinase that is implicated in the pathogen 
sensing signaling pathways in the immune system. 
IRAKI has been reported (Takeshita and Ishii 2008) to 
act as a triggered kinase in the detection of CpG DNA 
by Toll-like receptor 9 (TLR9). Besides, it was found 
that IRAKI plays a significant role as a signal 
amplifier that would culminate in the activation of 
NF-kB (O'Neill 2006). Although IRAKI was not found 
to involve in the inflammasome signaling pathway, 
IRAK4, a member from the same kinase family as 
IRAKI, was reported to assemble with MyD88 in the 
NLRP3 inflammasome activation (de Nardo and Latz 
2011 ). 

In the TLR signaling pathway, MyD88 adaptor protein 
is associated with TLRs via Toll-IL-1 receptor (TIR) 
domain interactions. The bound TLR-MyD88 complex 
then recruits IRAK4 by binding to its death domain. 
The activated IRAK4 further phosphorylates IRAKI, 
leading to the activation of downstream Traf-6, Tak-1, 
inhibitor of NF-kB (IKK), IRF7, and transcription factor 
NF-kB (O'Neill 2006; Blasius and Beutler 2010). 
Activated NF-kB will translocate into nucleus and 
regulate the expression of type I interferon production. 
The appropriate level of type I interferon is critical not 


only in the immune defense against pathogens (Axtell 
and Steinman 2008); its aberration will have enormous 
impacts on the autoimmunity (Ewald and Barton 2011). 

Structural insight of the proteins is important for the 
disclosure of the function. Besides, structural studies 
also provide clues on the interaction between ligands 
and their receptors (Wei et al. 2011). In the field of 
protein science, structural studies normally are carried 
out by electron microscopy (Nasertorabi et al. 2011), 
NMR techniques (Wang et al. 2011), and 
crystallography (Piotrowski et al. 2009). These 
approaches are complementary in the efforts to better 
understand the protein structure and function (Hunter 
et al. 2011). Despite the widespread use of physical 
approaches, computational approach has also been 
used increasingly due to its efficiency in the biological 
research. To date, computational analysis of the 
biological processes and entities has been applied in 
various fields, such as cellular pathways (Zhang et al 
2009), cytoskeleton studies (Dodson and 
Dimitrakopoulos 2010), protein structure (Kim 2004; 
Arakaki et al. 2004; Hinsen 2008), RNA (Levenkova et 
al. 2004; Jonikas et al. 2009; Muller 2005), viral studies 
(Zhou et al. 2005; Itakura et al. 2010), immunogenetics 
(Massari et al. 2012), signalling pathways (Xu and 
Wong 2008; Fujarewicz et al. 2007; Lai et al. 2009), 
cellular energies (Andronescu et al. 2010), oncology 
(Han 2010; Wang et al. 2007; Oh and Gao 2011; Wang 
et al. 2011), molecular docking (Axenopoulos et al. 
2011; Zhang et al. 2005; Hussain et al. 2010; Teixeira et 
al. 2011), protein folding (He et al. 2009; Roche et al. 
2011; Csaba and Zimmer 2010), carbohydrate structure 
(Ranzinger et al. 2011), and drug discovery (Sato et al. 
2012; Kushwaha and Shakya 2010). 

In this research, a computational approach has been 
adopted to analyze the domains of IRAKI protein. 
Both active site and substrate binding site of IRAKI 
were studied in detail. The best template for IRAKI 
was predicted in this study. The insights of this 
research may guide further research in the 
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immunotherapy studies of IRAKl-implicated 
pathways. 

Methods 


target protein and template protein was carried out 
using SCWRL (Krivov et al. 2009) and JACKAL 
modelling approach. Lastly, Peons server (Lundstrom 
et al. 2001) was used to validate the structure of IRAKI. 


The nucleotide and amino acid sequence of IRAKI 
were retrieved from National Center for 
Biotechnology Information (NCBI). We used PSI- 
BLAST (Altschul et al. 1997) for distant homology 
recognition. Cn3D version 4.3 was used to simulate 
the domains of the protein. FFAS03 method 
(Jaroszewski et al. 2005) has been used for the 
calculation of the confidence score for the structural 
domains of IRAKI. Sequence alignment between 


Results and Discussion 

Sequence alignment was performed on the active site 
and substrate binding site of IRAKI. Fig. 1 illustrates 
the alignment for active site, with the domain 
highlighted in yellow. All of the proteins shown in Fig. 
1 are kinases. Fig.2 illustrates the alignment for 
substrate binding site. 
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FIG. 1 SEQUENCE ALIGNMENT OF THE ACTIVE SITE 
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FIG. 2 SEQUENCE ALIGNMENT OF THE SUBSTRATE BINDING SITE. 
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The abbreviation of the molecules presented in Fig. 1 
and Fig. 2 is as such: 1F3M_C is Chain C, Crystal 
Structure of H u m a n Serine/Threonine Kinase PAK1; 
1L3R_E is hain E, Crystal Structure of A Transition 
State Mimic of The Catalytic Subunit of Camp- 
Dependent Protein Kinase; 1FGI_A is Chain A, Crystal 
Structure of The Tyrosine Kinase Domain of Fibroblast 
Growth Factor Receptor 1 In Complex With Su5402 
Inhibitor; 1 T4H_P) is Chain B, Crystal Structure of The 
Kinase Domain of Wnkl; gi 68125618 is protein kinase 
(Leishmania major strain Friedlin); gi 48428484 is 
Serine/threonine-protein kinase PINK1; gi 78045619 is 
protein kinase (Xanthomonas campestris pv. 
vesicatoria str. 85-10); gi 50057639 is Protein kinase, 
(Paramecium tetraurelia); gi 125291 is 
Serine/threonine-protein kinase HAL4/SAT4. 

From the sequence alignments, it was found that the 
extent of sequence identity is low, hence it is quite 
unlikely to identify a meaningful putative functional 
relationship, which is one of the objectives of the 
sequence alignment (Batzoglou 2005). The simulated 
domains are given in Fig. 3 (Ball and Stick model) and 
Fig. 4 (worms model), respectively. Homology 
relationship can be accurately recorded by the 
structural similarity, due to the fact that structural 
similarity is preserved beyond sequence similarity. 
Understanding of the three-dimensional structure of a 
protein may facilitate the understanding of the 
function of each domain (Vogel et al. 2004). 



FIG. 3 PROTEIN DOMAINS (BALL AND STICK MODEL) 



FIG. 4 PROTEIN DOMAINS (WORMS MODEL) 

FFAS03 method (Jaroszewski et al. 2005) has been 
used for the calculation of the confidence score for the 
structural domains of IRAKI. Table 1 lists the top 10 
best matched templates with IRAKI. It is apparent that 
2nru protein is the most probable protein template 
that matches the structure of IRAKI. 


TABLE 1 BEST MATCHED PROTEIN TEMPLATES FOR IRAKI 
(USING FFAS03) 


Confidence 

Score 

Targets 

Protein 

length 

(aa) 

Description 

-40.40 

2nru 

307 

Interleukin-1 receptor- 
associated kinase 4 

-40.00 

3hgk 

327 

Protein kinase 

-37.70 

3tl8 

349 

BRASSINOSTEROID 
INSENSITIVE 1- 

associated rece 

-30.80 

2qlu 

314 

Activin receptor type 
IIB 

-30.80 

3q4t 

322 

Activin receptor type- 
2A 

-30.40 

3p86 

309 

Serine/threonine- 
protein 
Kinase CTR1 

-30.10 

4f0f 

287 

Serine/threonine- 
protein 
Kinase roco4 

-29.80 

3g2f 

336 

Bone morphogenetic 
Protein receptor type- 
2 

-29.10 

3dtc 

271 

MAPK 

-28.00 

3s95 

310 

LIM domain kinase 1 


Based on the confidence score listed in Table 1, we 
identified 2nru as the potential protein template to 
model the structure of IRAKI. When there is a clear 
sequence relationship between a target protein and a 
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known structure, structural templates can be used to 
yield an accurate core model (Moult 2005). We 
modelled IRAKI against 2nru to find the identical 
residues. Such approach has been used in the 
experiment to study the glycosylase domain of DML 
proteins (Ponferrada-Marin et al. 2011). In this study, 
IRAKI was aligned with 2nru (aa: 15) starting at 
amino acid position 208. Fig. 5 demonstrates the 
identical residues of both strands, as highlighted in 
red. 

208 WPLCEISRGTHNFSEEL KIGEGGFGCVYRAVHRNTVYAVKRXXX 251 

IS FSFYELKNVTNNFDERPISVGGNKHGEGGFGWYKGYVNNTTVAVKK 61 

252 LKENADLEWTAVKQSFLTEVEQLSRFRHPNIVDFAGYCAQNGFYCLVYGF 301 

62 LAAMVDITTEELRQQFDQEIKVHAKCQHENLVELLGFSSDGDDLCLVYVY 111 

302 LPNGSLEDRLXXXHCQTQACPPLSUPQRLDILLGTARAIQFLHQDSPSLI 351 

112 MPNGSLLDRLSCL DGTPPLSWHMRCKIAQGAANGINFLHENH — HI 155 

352 HGDIRSSNVLLDERLTPKLGDFGXXXLARFSRFAGSSPSQSSMVARTQTV 401 

156 HRD IKSANILLDE AFTAKISDFGLARAS E KF AQTVMTS R I 195 

402 RGTLAYLPEEYIKTGRLAVDTDTFSFGWVLETLAGXXXQRAVKTHGART 451 

196 VGTTAYMAPE AL-RGE ITPKSD I YSFGWLLE IXTGLP AVDEHREPQLLL 244 

452 KYLKDLVEEEAEEAGVALRSTOSTLQAGLAADAWAAPIAHOIYKKHLDPX 501 

245 DIKEEIEDEEKT XEDYIDKK 264 

502 XXRPGPCPPELGLGLGQL AC CCL HRR AKRRP PHTQVYERLEKLOA 546 

265 M ND AD STS VE AMYS V AS QC L HE KKNKRP D I KKVQQL L QEMTA 306 

FIG. 5 PROTEIN STRUCTURE MODELLING 

To validate our results, we used Peons server to 
predict the structure of IRAKI. Table 2 lists the 
predicted matches. 


TABLE 2 BEST MATCHED PROTEIN TEMPLATES FOR IRAKI 
(USING PCONS) 


Rank 

Score 

Targets 

1 

3e-98 

2nru 

2 

3e-90 

lvjy 

3 

4e-89 

2b7a 

4 

4e-87 

2fb8 

5 

4e-86 

2eva 

6 

4e-86 

lyvj 

7 

le-80 

lsm2 

8 

9e-78 

2dq7 

9 

2e-77 

lk2p 


Comparing Table 2 with Table 1, 2nru is the best 
template yielded using FFAS03. Besides, it is also 
validated by Peons server. It has four chains (A, B, C, 
D) with a structure weight of 140266.73. The crystal 
structure of this protein in the phosphorylated form 
has been reported by Wang et al. 2006. Further 
evidence from the bound protein assembly suggests 
that both IRAKI and 2nru are assembled with MyD88 
in the similar way (Honda et al. 2004). Thus, 2nru 
represents the best matched protein template for 
IRAKI. 


Conclusion 

IRAKI is one of the important kinases in the pathogen 
sensing pathway. It is found to bind to the MyD88 
adaptor protein in the signaling cascade, which is 
leading to the activation of NF-kB transcription factor. 
The domain study of IRAKI is important for the 
disclosure of its functional role in the cytokine 
production and subsequent recruitment of 
lymphocytes to the sites of inflammation. In this study, 
we found that the structure of IRAKI is closely 
resembled to 2nru, which suggests that IRAKI and 
IRAK4 are highly homologous. Because both IRAKI 
and IRAK4 have domain similarity to MyD88, it is of 
interest to study other potential MyD88-bound 
proteins which are yet to discover. The domain 
analysis of IRAKI may provide insights into the 
therapeutic strategies that target MyD88 signaling 
pathway. 
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